February 5, 2008 


Constraints on the DGP Model from Recent Supernova 
Observations and Baryon Acoustic Oscillations 


Zong-Kuan Guo t, Zong-Hong Zhu °, J.S. Alcaniz 3, Yuan-Zhong Zhang *1 


ABSTRACT 


Although there is mounting observational evidence that the expansion of our 
universe is undergoing a late-time acceleration, the mechanism for this accelera- 
tion is yet unknown. In the so-called Dvali-Gabadadze-Porrati (DGP) model this 
phenomena is attributed to gravitational leakage into extra dimensions. In this 
work, we mainly focus our attention to the constraints on the model from the 
gold sample of type Ia supernovae (SNela), the first year data from the Supernova 
Legacy Survey (SNLS) and the baryon acoustic oscillation (BAO) peak found in 
the Sloan Digital Sky Survey (SDSS). At 99.73% confidence level, the combi- 
nation of the three databases provides Qm = 0.270+9:018 and Q, = 0.21673 313 
(hence a spatially closed universe with Q} = —0.350+9-980) which seems to be 
in contradiction with the most recent WMAP results indicating a flat universe. 
Based on this result, we also estimated the transition redshift (at which the uni- 
verse switches from deceleration to acceleration) to be 0.70 < zg-o < 1.01, at 20 
confidence level. 


Subject headings: cosmological parameters — cosmology: theory — distance 
scale — supernovae: general — galaxies:general 


1. Introduction 


Recent observations of type Ia supernovae (SNe Ia) suggest that the expansion of the 
universe is accelerating (Riess et al. 1998, Perlmutter et al. 1999, Tonry et al. 2003, Barris 
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et al. 2004, Knop et al. 2003, Riess et al. 2004). As is well known all usual types of 
matter with positive pressure generate attractive forces, which decelerate the expansion of 
the universe. Given this, a dark energy component with negative pressure was suggested to 
account for the invisible fuel that drives the current acceleration of the universe. There are a 
huge number of candidates for the dark energy component in the literature (see, e.g., Sahni 
and Starobinsky 2000; Peebles and Ratra 2003; Padmanabhan 2003; Lima 2004; Copeland 
et al. 2006 for recent reviews), such as a cosmological constant A (Carroll et al. 1992), an 
evolving scalar field (referred to by some as quintessence: Ratra and Peebles 1988; Caldwell 
et al. 1998; Weller and Albrech 2002; Guo, Ohta and Zhang 2005), the phantom energy, 
in which the sum of the pressure and energy density is negative (Caldwell 2002; Dabrowski 
et al. 2003; Wu and Yu 2005a), the quintom model (Feng, Wang and Zhang 2005; Guo 
et al. 2005; Zhao et al. 2005; Wu and Yu 2005b), the holographic dark energy (Li 2004; 
Gong 2004; Wang, Gong and Abdalla 2005; Myung 2005; Zhang and Wu 2005; Pavon and 
Zimdahl 2005; Chang, Wu and Zhang 2006), the Chaplygin gas (Kamenshchik et al. 2001; 
Bento et al. 2002; Dev, Alcaniz and Jain 2003; Silva and Bertolami 2003; Makler et al. 2003; 
Zhu 2004; Gong 2005; Zhang and Zhu 2006), and the Cardassion model (Freese and Lewis 
2002; Zhu and Fujimoto 2002, 2003; Godlowski, Szydlowski and Krawiec 2004; Amarzguioui, 
Elgaroy and Multamaki 2005; Koivisto, Kurki-Suonio and Ravndal 2005; Lazkoz and Leon 
2005; Szydlowski and Godlowski 2006). 


Another possible explanation for the accelerating expansion of the universe could be the 
infrared modification of gravity expected from extra dimensional physics, which would lead 
to a modification of the effective Friedmann equation at late times. An interesting model 
incorporating modification of gravitational laws at large distances was proposed by Dvali, 
Gabadadze and Porrati (2000), the so-called DGP model. It describes our four-dimensional 
world as a brane embedded into flat five-dimensional bulk. While ordinary matter fields are 
supposed to be localized on the brane gravity can propagate into the bulk. Unlike popular 
braneworld theories at the time, the extra dimension featured in this theory is astrophysically 
large and flat (for a recent review of the DGP phenomenology, see Lue 2005). A crucial 
ingredient of the model is the induced Einstein-Hilbert action on the brane. In this model, 
gravitational leakage into the bulk leads to the observed late-time accelerated expansion of 
the universe. Such a possible mechanism for cosmic acceleration has been tested in many of 
its observational predictions, ranging from local gravity (Lue 2003; Lue and Starkman 2003; 
Lue, Scoccimarro and Starkman 2004) to cosmological observations, such as SNe Ia (Deffayet, 
Dvali and Gabadadze 2002; Deffayet et al. 2002; Avelino and Martins 2002; Dabrowski et 
al. 2004; Alam and Sahni 2005; Maartens and Majerotto 2006), angular size of compact 
ratio sources (Alcaniz 2002), the age measurements of high redshift objects (Alcaniz, Jain 
and Dev 2002), the optical gravitational lensing surveys (Jain, Dev and Alcaniz 2002), the 
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large scale structures (Multamaki et al. 2003), and the X-ray gas mass fraction in galaxy 
clusters (Zhu and Alcaniz 2005; Alcaniz and Zhu 2005). 


This paper aims at placing new observational constrains on the DGP model by using 
the gold sample of 157 SNe Ia compiled by Riess et al. (2004), the 71 new SNe Ia released 
recently by the Supernova Legacy Survey (SNLS) (Astier et al. 2005), and the baryon 
acoustic oscillations detected in the large-scale correlation function of Sloan Digital Sky 
Survey (SDSS) luminous red galaxies (Eisenstein et al. 2005). It is shown that, if only 
the SNe Ia databases are used, 9,, and Qm are highly degenerated. However when we 
combine the baryon acoustic oscillations found by Eisenstein et al. (2005) from the SDSS 
data for analyzing, the degeneracy between Q,, and Qm is broken and the two parameters 
are accurately determined. 


We stuctured this paper as follows. Section 2 discusses the basic expressions of the DGP 
model. In Section 3, we present our analysis of the model using the updated SNe Ia data and 
the baryon acoustic oscillations found in the SDSS data. We end the paper by discussing its 
main results in Section 4. 


2. Basic expressions of the DGP model 


In the DGP model the modified Friedmann equation due to the presence of an infinite- 
volume extra dimension reads (Deffayet, Dvali and Gabadadze 2002; Deffayet et al. 2002) 


H? = Rè | (1+ 2) + ( 0, + On Fal =F) | (1) 


where H is the Hubble parameter (Hp is its current value), Q; and Qm represent the fractional 
contribution of curvature and of the matter (both baryonic and nonbaryonic), respectively, 
and Q,., the bulk-induced term, is defined as 


Qn, = 1/4r2 Hè. (2) 


In the above equatios, re is the crossover scale beyond which the gravitational force follows 
the 5-dimensional 1/r? behavior. Note that on short length scales r < re (at early times) 
the gravitational force follows the usual four-dimensional 1/r? behavior, i.e., the standard 
cosmological models are recovered. It has been shown that by setting the crossover scale re 
close to the horizon size, this extra contribution to the Friedmann equation leads to acceler- 
ation which can in principle explain the supernova data (Deffayet, Dvali and Gabadadze 
2002; Deffayet et al. 2002). From Eq. 1 we find that the normalization condition is 
given by Qk + (Vr + /Q,, + Qm)? = 1, while for a spatially flat scenario it reduces 
to Qr, = (1 0) Ay, 


s= 


The current value of the deceleration parameter, defined q = —aä/à?, takes the form 
(Zhu and Fujimoto 2003, 2004) 


3 V Qr 2 
= -Qm | L+ ——=— | — Q, + Vr Fm) . 


The transition redshift z4=0 at which the universe switches from deceleration to acceleration, 
can be expressed in the following analytic form (Zhu and Alcaniz 2005) 


Q 1/3 
an=) : (4) 


Note that from a phenomenological standpoint, the DGP model is a testable scenario 
with the same number of parameters as the ACDM scenario, contrasting with models of 
quintessence that have additional free parameters to be determined (Deffayet et al. 2002). 


3. Constraints from SNela and SDSS data 


In this section we analyze the DGP model by using two recently released supernova 
data sets, the Gold supernova data set (Riess et al. 2004) and the SNLS data set (Astier 
et al. 2005). We also use these data sets in conjunction with the recent discovery of the 
baryon acoustic oscillation peak in the SDSS (Eisenstein 2005) to place constrains on the 
cosmological parameters. 


Recently, Riess et al. (2004) compiled a large database of 170 previously reported SNe 
Ia and 16 new high redshift SNe Ia observed by the Hubble Space Telescope (HST). The 
total sample span a wide range of redshift (0.01 < z < 1.7). To reflect the difference in the 
quality of the spectroscopic and photometric record for individual supernovae, they divided 
the total sample into “high-confidence” (gold) and “likely but not certain” (silver) subsets. 
Here, we consider only the gold sample of 157 SNe Ia (for recent usages of the sample, see, 
e.g., Padmanabhan and Choudhury 2003; Nesseris and Perivolaropoulos 2004; Alcaniz 2004; 
Choudhury and Padmanabhan 2005; Gong 2005; Feng, Wang and Zhang 2005; Zhang and 
Wu 2005; Guo and Zhang 2005a,b; Cai, Gong and Wang 2006; Ichikawa and Takahashi 
2005). 


More recently, the SNLS collaboration released the first year data of its planned five- 
year Supernova Legacy Survey (Astier et al. 2005). An important aspect to be emphasized 
on the SNLS data is that they seem to be in a better agreement with WMAP results than 
the gold sample (see, e.g., Jassal, Bagla and Padmanabhan 2006). The two samples are 
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illustrated on a residual Hubble Diagram with respect to our best fit universe (Qm = 0.270, 
Q,, = 0.216) in Figure 1. 


It is well know that the acoustic peaks in the cosmic microwave background (CMB) 
anisotropy power spectrum can be used to determine the properties of the cosmic pertur- 
bations, to measure the contents and curvature of the universe, as well as many other cos- 
mological parameters (see, e.g., Spergel et al. 2003). Because the acoustic oscillations in 
the relativistic plasma of the early universe will also be imprinted on to the late-time power 
spectrum of the non-relativistic matter (Peebles and Yu 1970; Eisenstein and Hu 1998), the 
acoustic signatures in the large-scale clustering of galaxies yield additional tests for cosmol- 
ogy. In particular, the characteristic and reasonably sharp length scale measured at a wide 
range of redshifts provides distance-redshift relation, which is a geometric complement to 
the usual luminosity-distance from type Ia supernove (Eisenstein et al. 2005). Although the 
acoustic features in the matter correlations are weak and on large scales, Eisenstein et al. 
(2005) have successfully found the peaks using a large spectroscopic sample of luminous, red 
galaxies (LRGs) from the Sloan Digital Sky Sruvey (SDSS, York et al. 2000). This sample 
contains 46,748 galaxies covering 3816 square degrees out to a redshift of z=0.47. They 
found a parameter A, which is independent of dark energy models (Eisenstein et al. 2005). 
From their Eq. 2 and 4, we write it as follows, 


Ve | z 1 ? ( f dz J 
A = — |— sinn’ | yI — ; 5 
a (Ela) [Qe Mal | EE ©) 
where E(z) = A(z)/Ho, z1 = 0.35, A is measured to be A = 0.469 + 0.017, and the 
function sinn(z) is defined as sinn(x) = sin(x) for a closed universe, sinn(z) = sinh(«) for 
an open universe and sinn(«) = x for a flat universe. In our analysis, we will comebine these 
measuremets. 


In order to place limits on our Eq. (1), we perform a y?-statistics for the model param- 
eters (Qm, Qr.) and the Hubble constant Ho. Since we want to concentrate solely on the 
density parameters, we need to marginalize over the Hubble parameter Hp. However, Ho 
appears as a quadratic term in y? or, equivalently, appears as a separable Gaussian factor 
in the probability to be marginalized over. Thus marginalizing over Ho is equivalent to 
evaluating x? at its minimum with respect to Ho (Barris et al. 2004). Here, we marginalize 
over the Hubble parameter by using the analytical method of Wang et al. (2004). Fig- 
ure 2 shows the joint confidence contour at 68.3%, 95.4% and 99.7% confidence levels in 
the parametric space Qm — Qr, arising from the gold sample of SN Ia data and the SDSS 
baryon acoustic oscillations. The best-fit parameters for this analysis are Qm = 0.272 and 
Q,, = 0.211. Note that the best-fit value for Q,, leads to an estimate of the crossover scale 
re in terms of the Hubble radius Hj", i.e., re = 1.089Hj ‘. Compared to Figure 2 of Alcaniz 
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and Pires (2004), the model parameters are more tightly constrained by using the prior 
from the baryon oscillation results than by assuming a Gaussian prior on the matter density 
parameter, Qm = 0.27 + 0.04, as provided by WMAP team (Spergel et al. 2003). 


Figure 3 illustrates the allowed regions in the Om — Qr. plane by using the first year 
SNLS data in conjunction with the SDSS baryon acoustic oscillations (see also Fairbairn and 
Goobar (2005) for a similar analysis'). Our best-fit for this joint SNLS plus BAO analysis 
happens at Qm = 0.265 and Q,, = 0.216. The parameter space is considerably reduced 
relative to Figure 2 since the SNLS data set is more sensitive to the value of Q,, than the 
gold sample. 


In Figure 4 we show the joint confidence contours from the gold sample of SN Ia data 
and the first year SNLS data. In this case, the best-fit model happens for Qm = 0.31 and 
Q,.. = 0.23. We find that the degeneracies between these parameters are broken by combining 
these two data sets in the joint statistical analysis. With the prior from the SDSS baryon 
acoustic oscillations, our fits provide Qm = 0.270 and Q,, = 0.216. Compared to Figure 4, 
the allowed confidence regions are slightly reduced. 


Note that a closed universe is obtained at 30 confidence level in the above analyses, 
which confirms the previous results obtained using the SNe Ia and the X-ray mass fraction 
data of galaxy clusters (Zhu and Alcaniz 2005; Alcaniz and Zhu 2005). Although there is 
a range on the parameter plane being consistent with both the SNela and the SDSS data, 
and the resulting matter density Qm is reasonable, a closed universe is obtained at a 99% 
confidence level, which seems to be inconsistent with the result, Q, = —0.02+?03, found 
by the WMAP team (Bennett et al. 2003, Spergel et al. 2003) and Q, = 0 predicted by 
the simplest inflationary scenarios. Avelino and Martins (2002) analyzed the same model 
with the 92 SNe Ia from Riess et al. (1998) and Perlmutter et al. (1999). Assuming a flat 
universe, the authors obtained a low matter density and claimed the model was disfavorable. 
In additional to including new SN Ia data from, and combining the SDSS data, 
we relax the flat universe constraint in our analysis. We obtained a reasonable matter 
density, but a closed universe. This means that, in light of WMAP results — a nearly flat 
universe with Q, = —0.02+9:92 — the accelerating universe from gravitational leakage into an 
extra dimension seems not to be favored by the current observational data. Note also that the 
best fit values of Qm and Q,, lead to an estimate of the transition redshift zj-9 = 0.8619: 
which is larger than the one estimated from the gold sample, i.e., zą=0 = 0.46 + 0.13 (Riess 


et al. 2004). It means that acceleration in the DGP model happens earlier. Figure 5 shows 
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Table 1: Constrains on Qm, Qro, Te and Zg=0 


Gold Sample (sa 0.24003 1.02000 0.7812 
Gold+BAO 0.272+9-023 OS ay 1.089+9:970 0.840-11 
SNLS 0:23:33 O20" 1121979 Gore 
SNLS+BAO C2662 0,216+9013 1,076+9:935 0.874008 
Gold+SNLS CS ae: 2a 10a. Osi 


the deceleration parameter as a function of redshift z for our best-fit values in DGP model. 
For comparison, we have also plot the curve for the standard ACDM model. In Table 1 we 
summarize the main results of the paper. 


4. Conclusion and discussion 


Observations of SNe Ia indicate that the expansion of the universe is accelerating. What 
drives the acceleration, however, is still acompletely open question. From the observational 
viewpoint, it is of fundamental importance to differentiate between the two major possibili- 
ties, namely, the existence of new fields in high energy physics (dark energy) or modifications 
of gravitation theory on large scales. In this paper, we have focused our attention on one 
of the leading contender in modified-gravity explanation of acceleration, the so-called DGP 
model. We have analyzed the DGP model by using the gold SN Ia sample, the recent SNLS 
data and the SDSS baryon acoustic oscillations. Since SN Ia data are sensitive to the value 
of 2, while the baryon acoustic oscillations are sensitive to the value of Qm, the combina- 
tion of these data sets breaks the degeneracies between the model parameters and leads to 
strong constraints on them, as shown in Figures 2,3,4. The joint analysis strongly indicates 
a spatially closed universe, which was already obtained by fitting the combination of SN Ia 
data and the X-ray gas mass fraction in galaxy clusters (Zhu and Alcaniz 2005; Alcaniz and 
Zhu 2005). We also estimate the transition redshift z,29 > 0.70 at 20 confidence level. 


In summary, we have discussed the gravitational leakage into extra dimensions as an 
alternative mechanism for the late-time acceleration of the universe (and an alternative route 
to the dark energy problem). In agreement with other recent analysis, we have shown that 
a spatially closed DGP scenario with a crossover scale re ~ Ho" is largely favored by most 
of the current observational data. 
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Fig. 1.— The gold sample and the SNLS sample are shown in a residual Hubble diagram 
with respect to the DGP model with the best-fit parameters, (Om, Qr.) = (0.270, 0.216). 
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Fig. 2.— Probability contours at 68.3%, 95.4% and 99.7% confidence levels for Qm versus 
Q,,. in the DGP model from the gold sample of SNela data (solid contours), from the baryon 
acoustic oscillations found in the SDSS data (dotted lines) and from the combination of the 
two databases (coloured contours) — see the text for a detailed description of the method. The 
upper-left shaded region represents the “no-big-bang” region, the thick solid line represents 
the flat universe and accelerated models of the universe are above the the dashed line. The 
best fit happens at Qm = 0.272 and Q,, = 0.211. 
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Fig. 3.— Probability contours at 68.3%, 95.4% and 99.7% confidence levels for Qm versus Q., 
in the DGP model from the first year SNLS data (solid contours), from the baryon acoustic 
oscillations found in the SDSS data (dotted lines) and from the combination of the two 
databases (coloured contours). The upper-left shaded region represents the “no-big-bang” 
region, the thick solid line represents the flat universe and accelerated models of the universe 
are above the the dashed line. The best fit happens at Om = 0.265 and Q,,. = 0.216. 
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Fig. 4.— Probability contours at 68.3%, 95.4% and 99.7% confidence levels for Qm versus Q., 
in the DGP model from the combination of both the gold sample of SN Ia data and the first 
year SNLS data (solid contours), from the baryon acoustic oscillations found in the SDSS 
data (dotted lines) and from the conjunction of the three databases (coloured contours). The 
upper-left shaded region represents the “no-big-bang” region, the thick solid line represents 
the flat universe and accelerated models of the universe are above the the dashed line. The 
best fit happens at Qm = 0.270 and Q,, = 0.216. 
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Fig. 5.— The deceleration parameter as a function of redshift z for some best-fit values in 
DGP model and the standard ACDM. 


